Flux saturation number of superconducting rings 
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The distributions of electrical current and magnetic field in a thin-film superconductor ring is 
calculated by solving the London equation. The maximum amount of flux trapped by the hole, 
the fluxoid saturation number, is obtained by limiting the current density by the depairing current. 
The results are compare it with similar results derived for the bulk case of a long hollow cylinder 
[Nordborg & Vinokur, Phys. Rev. B 62, 12408 (2000)]. In the limit of small holes our result reduces 
to the Pearl solution for an isolated vortex in a thin film. For large hole radius, the ratio between 
saturation numbers in bulk and film superconductors is proportional to the square root of the hole 
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Use of thin film superconductors integrated in nanodc- 
vices requires precise knowledge of the film behavior in 
the presence of both external and self induced magnetic 
fields. Recent experiments have shown that properly de- 
signed arrays of dots and antidots can serve as effective 
traps for magnetic fluxpii2i£ It has also been shown that 
patterned superconducting films allow for the motion of 
magnetic vortices to be guided over the film area^&i^ 
opening up for a new field often called fluxonics. 

The models for flux trapping used in past years did not 
take into account the precise geometry of the patterned 
film samples. Those conventional models were based on 
approximations applicable cither to a single vortex in 
an infinite film, equivalent to having a hole of zero ra- 
dius, i.e., a Pearl vorte^, or to an infinitely long hollow 
cylinder . 10 ' 11 Both geometries are far from realistic for 
thin film devices. However several recent works consider- 
ing dynamics of vortices, current and field distributions in 
thin-film superconductor with finite hole size have been 
publishedi 12 ' 13 i 14 i 15 i 16 Yet, the flux saturation number, 
defined as the maximum number of flux quanta which 
can be trapped inside the hole, have so far not been cal- 
culated. Since the saturation number is a key concept 
for the vortex-antidot interaction and the overall effect 
of film patterning, we dedicate the present work to this 
calculation. 

In this paper, we examine the problem of trapped flux 
in a thin-film superconducting ring by solving the Lon- 
don equations. The ring geometry allows modelling of a 
realistic situation with an isolated anti-dot in a thin-film 
sample. We consider the case of zero external magnetic 
field, and show that the saturation number for a thin-film 
ring can differ significantly from that of a hollow cylin- 
der with the same radius. It is shown how the difference 
depends on the ratio between the hole radius and the 
thickness of the ring. 

Consider a ring where the outer radius, r%, is much 
larger than the inner radius, r±, see Fig. [1] It is as- 
sumed that the film thickness, d, is negligible compared 
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FIG. 1: The thin ring geometry. 



to both the hole radius and the effective penetration 
length A ff = X 2 /d, where A is the London penetration 
depth of the superconductor. We want to determine the 
amount of flux trapped in the ring in a remanent state 
with flux present due to some magnetic prehistory. It 
will be assumed that flux pinning elsewhere is absent. 

In the superconductor n < r < r<x, and —d/2 < z < 
d/2, the distributions of current density j and induction 
B are given by the London equation, which in terms of 
the vector potential A, where curl A = B, reads 



curl(A 2 /i j + A) = . 



(1) 



Due to the symmetry, the current density and vector 
potential have in cylindrical coordinates only azimuthal 
components, j = (0,7,0) and A = (0,^4,0), respectively. 
From the Eq. it then follows that 



HoX 2 j(r)+A(r) = §t, 



(2) 



where flux quantization implies that the constant 
the London fluxoid, is a free parameter restricted to 
an integer number of the flux quantum, = n<&o, 
where $o = h/2e. The second term on the left-hand- 
side represents the flux through the area of radius r; 



/J" B(r')2-Kr' dr' — A(r) 2nr, where B is the induction 
in the film plane, and is directed perpendicular to the 
plane. 

In the absence of an applied field, the vector poten- 
tial is only due to the current in the ring, and can be 
expressed as 
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(3) 



For the film geometry one can neglected variations of 
the current across the the sample thickness and average 
Eq. ([3]) over z. Introducing the sheet current J(r) = 

i( r )^ z i together with the dimensionless variables 

J(r) = OoA;y$/)J(r) and r = r/A cff the Eq. © be- 
comes, 



r~2 r2ir 



J(r') COS 9 



dOdr' 



n jo \/(f/f') 2 + 1 - 2(f/P) cos b 
= l/2vrf. (4) 

This Frcdholm integral equation of the second kind was 
solved numerically by converting it into a set of linear 
equations corresponding to discrete values of the radial 
coordinate, 
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o VWW + l-2(P/f,)cos# 



(5) 

dO. (6) 



Results of such calculations are presented in Fig. [51 The 
upper panel shows the field and current distribution for 
the case where r± = A e s , = 20A e ff. The field is mostly 
concentrated in the hole, where it has a pronounced peak 
at n, while a small peak from the return field is seen 
near ri. The shielding current flows predominantly in 
the vicinity of the hole, reaching a maximum at the edge. 

Plotted in the lower panel is the normalized cur- 
rent distribution for holes of different sizes ranging from 
fi/Xcs = 1 to 100, with all rings having V2/X e s = 1000. 
For comparison, the figure includes also the current dis- 
tribution around a single Pearl vortex in an infinite filing 



Jpeari(r) - [Si(f/2) - #i(f/2) - 2/tt] /8 . 



(7) 



Here Si is the first order Struve function and Ki the 
first order modified Bcsscl function of the second kind. 
For the small hole case, ri = A c ff, our numerical result 
is very close to the Pearl solution over the whole ring 
area except near the outer edge, where the sheet current 
has an upturn. The curves for the two larger hole radii 
also follow Eq. {7} for intermediate r, but are seen to 
have an additional upturn at the inner edge. The cur- 
rent enhancement near n is found to increase with the 
hole size, whereas near the outer edge the behavior is 
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FIG. 2: (Upper) Distribution of magnetic field H and current 
J in the superconducting ring. Inner radius of the ring n is 
equal to A c ff and outer radius r2 is equal to 20A c ff. (Lower) 
Current density distribution in superconducting thin rings 
with 3 different hole radii, n, and the same overall size ri. For 
comparison, the graph includes also the Pearl vortex current 
distribution, plotted as the dash-dotted line. 



marginally influenced by the hole. Similar results for su- 
perconducting rings were reported previously by Brandt 
and Clern^. 

Our main interest lies in finding the fluxoid saturation 
value as function of the size of the hole in a thin film. It 
follows from Eqs. @ and © that the magnitude of the 
sheet current depends linearly on $y. Thus, the max- 
imum <frf corresponds to having a current at the inner 
edge, J{ri), equal to the maximum current supported by 
the superconductor. We take this maximum value to be 
the depairing current^ 



Jdp 
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Thus, the fluxoid saturation number, n sat , satisfies, 



JdpC 



An) ' 



(8) 



(9) 



which gives 
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(10) 



Shown in Fig. [3] is the fluxoid saturation number as a 
function of the hole size. The following material parame- 
ters were here used, £ = 3 nm, A = 150 nm, d = 100 nm, 
corresponding to YBa2Cu30 x at low temperatures. The 
stepwise increase of n sat versus r\ seen in the main plot 
is essentially linear, but with an additional weak upward 
curvature. The inset shows that this behavior continues 
also as for larger holes, r\ > A. 

For small holes we found that the current near the 
inner edge behaves essentially according to Eq. ([7]) , which 
in this limit is well approximated by J(f) = l/2irr. Thus, 
for small holes the saturation number is given by, 



2ri 
3V3C 



ri < A e ff . 



(11) 




Hole size, r A 



FIG. 3: Trapped flux inside the hole of thin superconducting 
film (blue line) and inside the infinite cavity in bulk super- 
conductor (red line) depending on the size of the hole/cavity. 



For larger holes where r\ 3> A c ff, the current upturn near 
the edge prevents us from using the Pearl solution di- 
rectly. Instead, we make a conjecture that the sheet cur- 
rent distribution is described by 



J(f) 



(12) 



Firstly, this analytical form fits excellently when com- 
pared to the current enhancement calculated numerically. 
Secondly, this form has the same asymptotic behavior as 
the Pearl solution, J(r) = l/wr 2 for f > 1, consistent 
with the close agreement seen in Fig. [2] over a wide in- 
termediate range of r. Finally, the Eq. (fT2l) has the di- 
verging factor l/y/r — r\, commonly present in the edge 
behavior of thin superconductors in the Meissner stated 
As usual, the divergence is cut off at a distance A ff from 
the edge, giving for the ring geometry a maximum sheet 
current of J(fi) = l/(v / 2vrfl 3/ ' 2 ). It then follows from 
Eq. (fTTH) that for large holes one has, 
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Tl > A e ff . 



(13) 



It is of interest to compare these results with the flux- 
oid saturation number associated with a hole in a bulk 
superconductor. The case of a circular hole through an 
infinite superconductor was discussed in Ref. fTTL using 
that for this geometry with T\ being the hole radius, the 
induction is given by 



j{ r i) = jdp, the maximum number of trapped flux quanta 
becomes, 
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(15) 



The bulk saturation number is also plotted in Fig. [3l 
where one sees the same tendency as for the thin film 
case, namely that a larger hole size allows more flux to 
be trapped. For small r± the two curves are actually 
overlapping. This is because in Eq. (|2|), which is valid 
for any cylindrical geometry, the vector potential term 
becomes negligible for small r, giving /j,q\ 2 j(r)2irr = <£>/ 
as the asymptotic behavior independent of sample thick- 
ness. Thus, the linear hole size dependence in Eq. (jll[) 
holds also for the bulk case. 

As T\ increases the curves deviate, showing that thin 
film superconductors will not trap as much flux as a 
bulk sample of the same material. For ri/A 1, the 
Bessel functions in Eq. (TTS"]) can be simplified, and one 
gets asymptotically a quadratic hole size dependence, 
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n > A. 



(16) 



For large holes, the ratio in flux trapping capability is 
simply 
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(17) 



B(r) = n<S> CK (r/\) , r > r u 



(14) 



where 1/C = 7rriA[2A"i(ri/A) + (r 1 /X)K (r 1 /X)}, and 
Ao, K\ are the zeroth-order and first-order modified 
Bessel functions. The current density around the hole 
here obtained simply by, ^j{r) = B'(r), and with 



For typical patterned films with hole size^i r\ = 500 — 600 
nm, this ratio is close to 2. 

Note that when the current density approaches the de- 
pairing value, the superconducting order parameter will 
be reduced. In particular, this will occur most easily 
in the vicinity of the hole. This means that one should 
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consider the saturation numbers derived in this work as 
upper limiting values for the number of flux quanta pos- 
sibly trapped by the superconductor. 

In summary, we have derived the fluxoid saturation 
number for thin superconducting films containing a cir- 
cular hole of finite radius. For small holes, ri <C A e ff, 
the result is approximately equal to the flux saturation 
number for bulk superconductors with a cylindric cav- 
ity. For large holes, i\ 3> A c ff, the saturation number in 
thin films is less than in bulks by a factor y/rij2d. The 
present results, obtained by assuming the sample size to 
be much larger than the hole, should apply to supercon- 
ducting films patterned with a single antidot. Even more, 
knowing the fluxoid saturation number, one may predict 



the interaction between vortices and antidot arrays in 
various experimental configurations. 
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